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ABSTRACT 

The main goal of this work is the analysis of new approaches to the study of the 
properties of astronomical sites. In particular, satellite data measuring aerosols have 
recently been proposed as a useful technique for site characterization and searching 
for new sites to host future very large telescopes. Nevertheless, these data need to 
be critically considered and interpreted in accordance with the spatial resolution and 
spectroscopic channels used. In this paper we have explored and retrieved measure- 
ments from satellites with high spatial and temporal resolutions and concentrated on 
channels of astronomical interest. The selected datasets are OMI on board the NASA 
Aura satellite and MODIS on board the NASA Terra and Aqua satellites. A compar- 
ison of remote sensing and in situ techniques is discussed. As a result, we find that 
aerosol data provided by satellites up to now are not reliable enough for aerosol site 
characterization, and in situ data are required. 
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1 INTRODUCTION 

Most aerosols reaching the Canary Islands are either ma- 
rine, CINa, cryogenic emissions or of African (Sahara and 
Sahel) origin The latter (clays, quartzes, feldspars and cal- 
cites), because of their size, can reduce visibility in the opti- 
cal wavelength range and can therefore affect astronomical 
observations. 

Furthermore, aerosols cause radiative forcing, oceano- 
graphic deposits by winds (together with Fe and Al), 
nutrients and minerals for algae (a coastal increase in 
chlorophyll — phytoplankton biomass), sanitary effects, etc. 
Aerosols also play an important role in astronomical site 
conditions, producing more stable condensation nuclei, de- 
laying precipitation and causing the extinction, absorption, 
diffusion and reflection of extraterrestrial radiation. 

Most of the airmass flux component reaching the Ca- 
narian archipelago comes from the North Atlantic Ocean 
and consists of sea aerosols, which absorbs chloride in the 
UV. African dust intrusions affect the western and eastern 
Canary Islands differently. Moreover, the presence of a stable 
inversion layer and the pronounced orography of the western 
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islands (La Palma and Tenerife) produce different mass flux 
patterns in the lower (mixing) layers closer to the sea (up 
to 800 m) and in the median-upper (or free) troposphere 
layer (above the thermal inversion layer, i.e. above 1500 m), 
causing a seasonally dependent vertical drainage of airborne 
particles. 

There are remarkable differences between summer in- 
trusions, which can rise to the peaks of the mountains (high 
level gloom), at 2400 m, and those of winter, more frequent 
in the lower troposphere (anticyclonic gloom). 

Anticyclonic gloom is associated with strong, stationary 
anticyclonic conditions forced by dust accumulation between 
the soil and the inversion layer. It can favour the decrease in 
height of the inversion layer and the formation of clouds that 
do not easily precipate as rain, hence persisting for a longer 
time and providing a more stable sea of clouds. The dust is 
trapped by the sea of clouds (sea of dust) and is prevented 
from reaching the topmost level in the islands (above 1500 
m). 

The aerosol index provided by the TOMS (Total 
Ozone Mapping Spectrometer) is one of the most widely 
accepted products for detecting the daily aerosol content. 
TOMS Level 3 data are gridded in squares of 1° x 1.25° 
(latitude and longitude respectively) and are available 
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online at http://toms.gsfc.nasa.gov The spatial coverage 
is global 90°S-90°N and the temporal resolution is daily. 
Moreover, several techniques have been developed in situ 
to characterize the presence of dust locally at the Canarian 
observatories. In particular, a parameter related to sky 
transparency, the atmospheric extinction coefficient in the 
V (551 nm) and r (625 nm) bands has been measured at 
the Observatorio del Roque de los Muchachos (ORM) on 
La Palma since 1984 by the Carslberg Automatic Meridian 
Circle (CAMC). The archive is in the public domain at 
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http://www.ast.cam.ac.uk/~dwe/SRF/camc _extinction.htmll 
and provides a good temporal comparison with the values 
retrieved with remote sensing techniques from the Total 
Ozone Mapping Spectrometer on board the Nimbus7 
satellite and from other probes (Aura/OMI, Terra/MODIS 
and Aqua/MODIS, MSGl(Met8)/SEVIRI and EN- 
VIS AT/SCI AM ACHY). Our main aim is to overlap the 
geographical area of the ORM with the satellite data; for 
this reason we have used Level 2 data. Level data are 
the raw data from the satellite; Level 1 data are calibrated 
and geolocated, keeping the original sampling pattern; 
the Level 2 data used in this paper are converted into 
geophysical parameters but still with the original sampling 
pattern; finally the Level 3 data are resampled, aver- 
aged over space, and interpolated/averaged over time (from 




http:/ /people. cs.uchicago.edu/~yongzh/papers/CM_InXg_Scale. 



On examination, the Level 2 data have the same spa- 
tial resolution as the Instantaneous Field of View (IFOV) 
of the satellite. Through a software procedure, it is possible 
to create files containing information on geophysical vari- 
ables (data describing the solid earth, marine, atmosphere, 
etc., properties over a particular geographical area) and field 
values such as seconds, latitude, longitude, reflectivity in dif- 
ferent channels, the ozone column, the aerosol index, aerosol 
optical depth, cloud land and ocean fraction, SO2 and radi- 
ance. From remote sensing and in situ data it is possible to 
trace back to the cloud coverage and climatic trend. 



The purpose of this study is the analysis of new ap- 
proaches to the study of the aerosol content above astro- 
nomical sites. Our objective is to calibrate the extinction 
values in the V band (550 nm) (more details in section 3.1) 
with remote sensing data retrieved from satellite platforms. 



The paper is organized as follows: Section 2 describes 
the problem and background; Sections 3 and 4 concern the 
comparison of in situ atmospheric extinction data with the 
aerosol index provided by TOMS (Level 3 data); Sections 
5 and 6 deal with the analysis of Level 2 data from other 
satellites and their validity for site characterization by com- 
paring the satellite results with in situ measurements; and 
the summary and outlook are given in Section 7. 



Two appendices have been included with complemen- 
tary information. Appendix I describes the format of the 
satellite data, indicating the official websites for data ac- 
cess, and Appendix II includes a list of acronyms to aid in 
following the terminology used in the paper. 



Figure 1. Trade wind behaviour in the low mixing mar- 
itime layer (MML) and in the tropospheric layer (TL). From 
http: / / cip.es/personales/oaa/nubes/nubes.htm 



2 METEOROLOGICAL AND GEOPHYSICAL 
SCENARIOS 

2.1 The trade wind inversion as a determining 
factor of aerosol distribution 

Site testing campaigns are at present performed within the 
classical scheme of optical seeing properties, meteorological 
-paf^iSifefeis^dsi? (transparency, sky darkness, cloudiness, etc. 
New concepts related to geophysical properties (seismicity 
and microseismicity) , local climate variability, atmospheric 
conditions related to the optical turbulence (tropospheric 
and ground wind regimes) and aerosol presence have re- 
cently been introduced in the era of selecting the best sites 
for hosting a new generation of extremely large telescopes 
(which feature a filled aperture collector larger than 40 m, 
and which are considered worldwide as one of the highest pri- 
orities in ground-based astronomy), telescope and dome de- 
signs, and for feasibility studies of adaptive optics (Munoz- 
Tunon 2002; Munoz-Tunon et al. 2004; Varela et al. 2002; 
Varela et al. 2006; Munoz-Tunon et al. 2007). 

The Canarian Observatories are among the top sites 
for astronomical observations and have been monitored and 
characterized over several decades (Vernin et al. 2002). The 
trade wind scenario and the cold oceanic stream, in combi- 
nation with the local orography, play an important role in 
the retention of low cloud layers well below the summits to 
the windward (north) side of the islands, above which the 
air is dry and stable (the cloud layers also trap a great deal 
of light pollution and aerosols from the lower troposphere) . 

The trade winds and the thermal inversion layer be- 
tween 1000 and 1500 m (shown in Fig. Q~} have been the 
object of many studies over the last 50 years, either indi- 
rectly from observations of the stratocumulus layer (known 
locally as the "sea of clouds" ) , which forms at the condensa- 
tion level (Font-Tullot 1956) or from radiosondes (Huetz-de- 
Lemps 1969). For much of the year, the trade wind inversion 
layer separates two very different airmasses: the maritime 
mixing layer (MML) and the free troposphere layer (TL) 
(Torres et al. 2001). The presence of the inversion layer is 
crucial in the airmass flows reaching the islands and in the 
downflow at high elevations that they undergo. Saharan dust 
invasions mostly affect the eastern islands, but can occasion- 
ally reach the islands of Tenerife and La Palma; normally, 
however, they do not reach the level of the Observatories 
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(the ORM is 2396 m and the OT 2390 m above mean sea 
level) . 



2.2 Airmass intrusions in the Canary Islands 

Airmasses are classified into three types according to their 
origin and permanence over continental landmasses and sea: 
Atlantic (marine, ON a), European (anthropogenic emis- 
sions from sulphates and carbon) and African (dust, mineral 
aerosols). The most frequent are chlorides, i.e. clean salts of 
oceanic origin that do not affect astronomical observations, 
and that reach or exceed half the total contribution of mass 
flux over the summits. European airmasses are almost al- 
ways of anthropogenic origin (sulphates and carbon) and 
are of scant importance (between 0.8 and 7.2%; see Romero 
2003). 

A study of the origins and distribution of air flows 
performed by Torres et al. (2003) distinguished the ori- 
gins of airmasses that affect the lower (MML) and upper 
(TL) layers, the latter being at the height of the astro- 
nomical observatories. Daily isentropic retrotrajectories (at 
00:00 and 12:00 GMT) were used during the AEROCE (At- 
mosphere/Ocean Chemisty Experiment) at Izaha (2367 m 
above MSL, 1986-97) and the Punta del Hidalgo Lighthouse 
(sea level, 1988-97), each airmass being assigned an origin: 
North American, North Atlantic, Subtropical Atlantic, Eu- 
rope, Local and Africa, the summer and winter periods being 
separated. 

The airmass provenance in the lower (mixing) layer 
close the sea (MML) is Northern Atlantic (the mean is 
59.6%), European (the mean is 19%) and African (0% 
in summer and 23% in winter), whereas the troposphere 
layer (TL) is mainly Northern Atlantic (44.2%) and African 
(17.4%), with a minimum in April (5.3%) and a maximum 
in August (34.5%) (Torres et al. 2003). 

African summer intrusions are therefore almost absent 
in the MML and more intense in the TL because of the 
daily thermal convection reaching the higher atmospheric 
layers. In winter intrusions into the TL are less frequent. 
The airmass is carried horizontally by the prevailing wind 
and is affected by a process of separation, the larger particles 
(> 10 fj,m) leaving sediments at ground level over a short 
time and the smaller ones being carried across the Atlantic 
Ocean to distances of hundreds or thousands of kilometres 
from their place of origin. This sand creates a large feature 
(plume) that is visible in satellite images (see, for example, 
the EP/TOMS or Aura/OMI websites) and extends from 
the African coast across a band about 20° in latitude. During 
winter, the prevailing wind carries the dust from the south 
of Canary islands through an average of 10° in latitude to 
the Cape Verde Islands. Considering that the dust plumes 
reach an extension of about 2000 km, during the winter they 
rarely reach La Palma. Instead, in summer, winds above 
4 km in altitude can take these particles as far north as 
30° in latitude. In these conditions the dust plume over the 
ORM is composed mainly of small quartz particles in the 
range 0,5-10 /jm and the biggest particles precipitate (dry 
deposition). Typical dust storms take 3-8 days to disperse 
and deposit 1-2.4 million tonnes/year. Finally, the clouds 
of aerosols dissipate through advective processes or through 
rain (sumid deposition). 



3 DATABASE FOR THE ANALYSIS: IN SITU 
EXTINCTION COEFFICIENT VALUES AND 
AEROSOL INDEXES PROVIDED BY THE 
TOMS AND OMI SPECTROGRAPHS ON 
BOARD SATELLITES 

3.1 Atmospheric extinction coefficient measured 
above the ORM (La Palma) 

Atmospheric extinction is the astronomical parameter that 
determines transparency of the sky. Extinction is associ- 
ated with the absorption/scattering of incoming photons by 
the Earth's atmosphere and is characterized by the extinc- 
tion coefficient, K. Sources of sky transparency degradation 
are clouds (water vapour) and aerosols (dust particles in- 
cluded). This coefficient is wavelength-dependent and can 
be determined by making observations of a star at different 
airmasses. For details of the astronomical technique for de- 
riving the extinction coefficient values we refer the reader to 
King (1985). 

Long baseline extinction values for the 
ORM have been measured continously at the 
Carlsberg Auto matic Meridian Circle (C AMC, 
http:/ /www. ast.cam.ac.uk/dwe/SRF/camc_extinction. html) 



in the V band (550 nm) and more recently in the Sloan 
r' band (625 nm). To our knowledge, this is the largest 
available homogeneous database for an observing site. 

Extinction values and their stability throughout the 
night are essential for determining the accuracy of astro- 
nomical measurements. As nights with low and constant ex- 
tinction are classified as photometric, this parameter is con- 
sidered among those relevant for characterizing an observing 
site. 

On photometric dust-free nights the median of the ex- 
tinction is 0.19 mag/airmass at 480 nm, 0.09 mag/airmass 
at 625 nm and 0.05 mag/airmass at 767 nm. The extinction 
coefficients reveal that on clear days (denoted coronal-pure) 
the extinction values at 680 nm are below about 0.07-0.09 
mag/airmass, while on dusty days (diffuse-absorbing) they 
are always higher (Jimenez & Gonzalez Jorge 1998). The 
threshold that identifies the presence of dust is at 0.153 
mag/airmass (Guerrero et al. 1998). 

At the ORM, the extinction in V is less than 0.2 
mag/airmass on 88% of the nights, and extinction in ex- 
cess of 0.5 mag/airmass occurs on less than 1% of nights. A 
statistical seasonal difference is also detected (see Guerrero 
et al. 1998). 

Figure [2] shows the cumulative frequency of extinction 
over the ORM during winter (October-May) and summer 
(June-September). In summer, 75% of the nights are free 
of dust, while at other times of the year over 90% of the 
nights are dust-free (Guerrero et al. 1998) . These results are 
consistent with those provided by Torres et al. 2003. 



3.2 Characterization of aerosols in the near-UV 

Data from Total Ozone Monitoring Spectrometer (TOMS) 
and from the Ozone Monitoring Instrument (OMI) are anal- 
ysed to detect absorbing and non-absorbing aerosols at ul- 
traviolet (UV) wavelengths. A detailed description of these 
instruments and the products are given at the official Web- 
site (http://toms.gsfc.nasa.gov/). Absorbing aerosols in- 
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Figure 2. Frequency of extinction over the ORM during winter (top) and summer (centre). In both cases the modal value is 0.11 
mag/airmass. Their corresponding cumulative frequencies are also shown (bottom). The vertical line indicates the extinction coefficient 
limit for dusty nights, KV 0.153 mag/airmass (Guerrero et al. 1998). 



elude smoke deriving from biomass burning, industrial activ- 
ity, mineral dust, volcanic aerosols and soot. Non-absorbing 
aerosols are mostly sulphates. The UV spectral contrast is 
useful for retrieving values over land and ocean because of its 
low reflectivity in the spectrometer range. Backscattered ra- 
diation at A 340, 360 and 380 nm is caused mainly by molec- 
ular Rayleigh scattering, terrestrial reflection and diffusion 
by aerosols and clouds (through Mie scattering). Quantita- 
tively, aerosol detection from TOMS and OMI is given by: 

AjV A = -10o{log 10 [(^i) l-log 10 [fc) 1}, 

I- L\J360/meaaJ L\ 1 360 / calc J > 

where I mea s is the backscattered radiance at the wavelenghts 
measured by TOMS and OMI (with Mie and Rayleigh scat- 
tering, and absorption) and J ca i c is the radiance calculated 
from the model molecular atmosphere with Rayleigh scatter- 
ers. depends strongly on the absorbing optical thickness 
of the Mie scatterers. AN\ is also called the aerosol index 
(AI). AI > indicates the presence of absorbing aerosols 
and clouds (AI±0.2) and negative AI values indicate non- 
absorbing aerosols (Herman et al. 1997; Torres et al. 1998). 
Significant absorption has been set at AI > 0.7 by Siher et 
al. 2004 and at AI > 0.6 in this work (see explanation in 
next section and in Fig. 1 10 p . 

The next section (Results I) concerns the first compar- 
ison of the in situ atmospheric extinction coefficient with 
the aerosol index provided by TOMS using Level 3 data. In 
Section 5 we present recent instruments on board satellites 
suitable for aerosol content monitoring and in Section 6 (Re- 
sults II) we compare the atmospheric extinction coefficient 
with the aerosol index and aerosol optical depth provided 
by OMI and MODIS respectively, using Level 2 data. 



4 RESULTS I: COMPARISON OF AI 

PROVIDED BY TOMS WITH KV FROM 
THE CAMC 

Atmospheric extinction is related to the internationally rec- 
ognized term aerosol optical depth (AOD) (or thickness) 
and to the aerosol index. In this section we shall analyse 
the first results of comparing the atmospheric extinction 
coefficient with the aerosol index provided by TOMS (To- 
tal Ozone Mapping Spectrograph) on board NASAs Earth 
Probe satellite. We use Level 3 aerosol data, which are avail- 
able at: 

ftp://toms.gsfc.nasa.gOv/pub/eptoms/data/aerosol/l 



In a previous paper (Varela et al. 2004a) we presented 
the aerosol index from the sector corresponding to the Teide 
Observatory (Tenerife) against the atmospheric extinction 
coefficient recorded at Roque de los Muchachos Observatory 
(La Palma). The Teide Observatory (OT) is situated 2390 m 
above sea level in Izaha and the geographical coordinates are 
16° 30' 35" West and 28° 18' 00" North. The Observatorio 
del Roque de los Muchachos is situated on the edge of the 
Caldera de Taburiente National Park, 2396 m above sea level 
and the geographical coordinates are 17° 52' 34" West and 
28° 45' 34" North. The two observatories are about 133 km 
apart. We compared the AI from EP/TOMS data centred 
at the OT and at the ORM sectors (see figure 5 of Varela et 
al. 2004b). 

The consistency of AI from both boxes centred on the 
OT and on the ORM shown in Fig. points to a similar 
tropospheric aerosol distribution at both observatories. 

In the next section (Results II) we show the results of 
comparing the atmospheric extinction coefficient with the 
aerosol index provided by OMI on board NASA's Aura satel- 
lite and with the aerosol optical depth provided by MODIS 
on board NASA's Terra and Aqua satellites. 

To demonstrate why there is not necessarily any corre- 
lation between the AI and atmospheric extinction we have 
selected an intense invasion of African dust over the Canary 
Islands that occurred on 2000 February 26. Figure [4] shows 
the plume of dust as recorded by TOMS on board the Earth 
Probe. We have used AI values obtained from TOMS for the 
Canarian Observatories and atmospheric extinction values 
provided by the CAMC at the ORM during the night be- 
fore and after this episode. AI reaches its maximum (2.4 at 
the OT and 2.7 at the ORM) on February 27, precisely on 
the day when the plume reached Tenerife. Nevertheless, the 
CAMC measures an extinction value of less than 0.2 mag, 
with a high number of photometric hours. The reason is that 
when the plume arrived at La Palma it did not reach the 
level of the Observatory. 

In Fig. [S]we represent the atmospheric extinction in the 
V band (KV) provided by the Carlsberg Meridian Circle 
against the aerosol index provided by EP/TOMS from 1996 
to 2004. 

We have classified four quadrants using KV=0.15 
mag/airmass as the threshold for dusty nights and AI=0 
as the threshold of the presence of absorbing aerosols. 

Figure[5]shows a large number of TOMS data indicating 
the presence of absorbing aerosols coincident with CAMC 
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Figure 3. Aerosol index provided by EP/TOMS from 1996 to 2004 at the OT (red triangles) and at the ORM (blue crosses). These 
profiles indicate similar aerosol tropospheric distributions at both observatories. 



values that show low or no atmospheric extinction (bottom 
right in Fig. [5} . This result is due to the presence of a layer 
of dust below the Observatory level. This dust layer is high 
and/or thick enough to be detected by the TOMS. This 
condition appears in 17% of all cases. 

The case of agreement between large extinction coeffi- 
cient and large Al (top right in Fig.[5j is associated with dust 
presence in the upper troposphere layer (TL). We have ver- 
ified that this case occurs in only 11% of all measurements, 
with 58% of these corresponding to the summer months 
(June-September), when the warmest surface winds sweep 
dust from the African continent and rise towards the up- 
per layers by convective processes (Romero & Cuevas 2002; 
Torres et al. 2003). 

The opposite case, i.e. low extinction coefficient and low 
aerosol index (bottom left in Pig. [5J happens in 59% of cases. 

In the top left quadrant in Fig. [5] we find low Al val- 
ues of aerosol index but large KV (13% of the data points). 
A possible explanation is given by Romero & Cuevas 2002, 
who argue that the cause could be local and concentrated 
dust in a small area intercepting the light from a star (and 
therefore measured by CAMC) but that is not representa- 
tive of the average values from a 1° x 1.25° region provided 
by TOMS; or is due to the presence of high fine dust lay- 
ers above the Observatory that are not detected by TOMS. 
Here we propose cirrus or other clouds as the explanation; 
in fact, 78% of the points located in the top-left quadrant 
correspond to the winter months when the possibility of the 
presence of medium-high clouds is greater. 

In the four cases (four quadrants), the correlation coef- 
ficient and the square of the coefficient of the Pearson cor- 




Figure 4. Dust plume obtained from TOMS data 
over the western Sahara Desert and extending 
over the Atlantic Ocean and Canary islands. From 
http: / / toms . gsfc. nasa.gov / aerosls / africa/canary. html 



relation is smaller than 0.1, so there is no correlation at 
all between both parameters. The maximum correlation is 
found in the upper-right quadrant (mathematically, this is 
the first quadrant), i.e. Al positive and extinction coeffi- 
cient larger than 0.15. When this interval is narrowed to Al 
larger than 0.7 and KV larger than 0.2 mag/airmass (only 
4% of cases), this correlation increases slightly. Only when 
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Figure 5. Aerosol index provided by TOMS against the atmospheric extinction coefficient in V (integrated values) from the CAMC 
above the ORM. The red horizontal line is a threshold indicating the presence of dust in the atmosphere (Guerrero et al. 1998). The 
blue vertical dashed line is the threshold line indicating the presence of absorbing aerosols in the atmosphere. The top-right quadrant 
corresponds to points seen as dusty (from CAMC KV > 0.15) and of high absorbing aerosol presence (Al > 0.7 has been used by Siher 
et al. 2004). 



the summer period is considered under these last conditions 
can the correlation reach 0.55. 

This last result coincides with Siher et al. 2004 (see 
Fig. 2 in this paper), correlating Al TOMS data recorded 
by Nimbus7 and the atmospheric extinction coefficient at 
the ORM during summertime dusty events (AI>0.7 and 
KV>0.2). The correlation found is 0.76. It is important, 
however, to emphasize that this relatively high (0.76) cor- 
relation occurs for a subsample (4%) of the total cases that 
occur when AI>0.7 and KV>0.2. As result of the lack of sta- 
tistical meaning in this result, the yearly average extinction 
derived from TOMS (figure 6 in Siher et al. 2004) does not 
reflect the real values measured in situ at the Observatory. 

This correlation occurs in summertime, when the 
warmest winds can drive the dust to the level of the Ob- 
servatory. 

The non-correlation between Al satellite remote sensing 
data and in situ extinction data is also reported by Romero 
& Cuevas 2002 (see figure 3 in this paper), who compare 
the aerosol index (Al) provided by TOMS against the opti- 
cal aerosol depth (AOD) obtained with a multifilter rotating 
shadow band radiometer (MFRSR) at Izana Atmospheric 
Observatory (OAI) (close to the Teide Observatory, 133 km 
the ORM), an institution belonging to the Instituto Na- 
tional de Meteorologfa (INM), now called Agenda Estatal 
de Meteorologfa (AEMET). The MFRSR measures global 
and diffuse radiation in six narrow band channels between 
414 nm and 936 nm (with bandwidths between 10 and 12 
nm). 

Again, in Romero & Cuevas 2002 a low correlation coef- 
ficient is attained when the AOD is larger than 0.1 at 414.2 
nm and the Al is larger than 0.5 — coinciding with dust in- 
vasions taking place owing to convective processes from the 
African continent towards high levels — and will then be de- 
tected by both TOMS and the MFRSR. Otherwise, in dust- 



free conditions, there is no correlation between the in situ 
measurements (AOD) and AT 

The prevailing trade winds, the inversion layer presence 
and the abrupt orography of the western Canary islands (La 
Palma and Tenerife) play an important role in the downflow 
of dust at high elevations. 

The low spatial resolution of the EP/TOMS for astro- 
nomical site evaluation explains the absence of correlation; 
this spatial resolution cannot distinguish local effects since 
it averages over a surface equivalent to an entire island (a 
1° x 1.25° region given for TOMS Level 3 data made avail- 
able in GSFC TOMS archives, i.e. Ill km x 139 km); nei- 
ther does it distinguish the vertical dust drainage. This ef- 
fect is known as anticyclonic gloom, in many cases, even in 
such dust storms episodes, the airborne dust particles do not 
necessarily reach the level of the Observatory (more expla- 
nations in Section 2 of this paper). 

It is possible to retrieve high resolution data from 
TOMS instrument also with the spatial resolution of the 
IFOV (35 km x 35 km) or TOMS Level 2 data (see Bertolin 
2007) in order to show that the correlation between CAMC 
KV and EP/TOMS Al for the ORM site on La Palma shows 
a little improvement when data are analysed at higher spa- 
tial resolution. 

The next step in our study is to explore the use of 
other instruments on board different satellites that oper- 
ate in bands of astronomical interest (the visible and NIR) 
with higher spatial resolution than TOMS and with long- 
term databases (longer than a few years) . These instruments 
will provide updated high resolution images for astronomi- 
cal site evaluation (to date, only TOMS data have been used 
for aerosol content characterization above an astronomical 
observatory) . 
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5 RECENT INSTRUMENTS ON BOARD 
SATELLITES FOR AEROSOL CONTENT 
MONITORING FOR ASTRONOMICAL SITE 
EVALUATION 

5.1 Updated images 

We have explored the use of other detectors on board differ- 
ent satellites that operate in bands of astronomical interest 
(the visible and NIR) with higher spatial resolution than 
TOMS. 

The new generation of satellites Terra&Aqua/MODIS, 
Aura/OMI, MSG1/SEVIRJ, etc. — acronyms described in 
Appendix II — provide better-resolution images that can be 
used to demonstrate the existence or otherwise of any cor- 
relation between the presence of aerosols and atmospheric 
extinction. Similarly to the case mentioned above (February 
2000), we examine 2007 March 10, when a thick plume of 
dust blew over the Canary islands from the west coast of 
Africa. 

The Moderate Resolution Imaging Spectroradiometer 
(MODIS) on NASA's Aqua satellite took the picture shown 
in Fig. [6] for this day (with much better resolution than 
TOMS, see Table 1). The MODIS Aerosol (MOD04_L2) 
product monitors the ambient aerosol optical thickness over 
the oceans globally and over a portion of the continents and 
contains data that have a spatial resolution (pixel size) of 10 
km x 10 km (at nadir). More information about grids and 
granule coverage is available in the official MODIS website 
(http://modis.gsfc.nasa.gov/). We can see that the eastern 
islands are more affected by the dust plume than the west- 
ern ones. Also, the abrupt orography and the inversion layer 
play an important role in retaining the dust below the sum- 
mits (where the observatories are located) of the highest 
western islands (Tenerife and La Palma). This effect is the 
above-mentioned anticyclonic gloom. 

This situation has been verified by data provided by the 
CAMC that indicate no atmospheric extinction at all: 

• Extinction Coefficient in r'= 0.083 (good quality dust- 
free night) i.e. KV = 0.12 mag/airmass. 

• Number of hours of photometric data taken=10.22 

• Number of hours of non-photometric data taken=0.00 

In the next section we summarize data retrieved from 
the new generation of satellites. Despite the much better 
spatial and spectral resolution of the recent satellite aerosol 
measurements (AI and AOD), they are not at the moment 
sufficient for the aerosol content monitoring of an astronom- 
ical observatory. In situ data are also required, in particular 
at those astronomical sites with abrupt orography (ORM, 
Mauna Kea or San Pedro Martir). Spatial resolution of the 
order of the observatory area is needed. 

5.2 Updated data 

In order benefit from satellite data for local site characteriza- 
tion, we have gathered and studied NASA and ESA satellite 
data planned to retrieve information about aerosol, clouds, 
ozone and other trace gases (N2, O2, H2O, CO2, CH4) that 
are found in the terrestrial atmosphere. In this paper we 
have centred our analysis on the aerosol content. 

An overview of selected satellites, parameters and sam- 
pled periods is given in Fig. [7] which includes other pa- 



rameters (ozone, cloud condensei nuclei (CCN) and cloud 
fraction) to be analysed in a future paper. 

To ensure that the retrieved remote sensing data fields 
from different satellites (such as aerosol values or geolocation 
parameters) are precisely over the ORM site coordinates to 
compare with the atmospheric extinction by CAMC, we de- 
cided to work with Level 2 data, which have a projected 
effective pixel size given by the instantaneous field of view 
(IFOV). We have also selected the longer term database 
(to retrieve more data for the KV comparison). The satel- 
lites/instruments and parameters used in this work are sum- 
marized in Table 1. 

The previous parameters selected and retrieved are the 
aerosol index (AI) provided by OMI (Ozone Monitoring In- 
strument) on board Aura — with visible and ultraviolet chan- 
nels and with a spatial resolution from 13 km x 24 km to 24 
km x 48 km — and the aerosol optical depth (AOD) provided 
by MODIS (Moderate Resolution Imaging Spectroradiome- 
ter) on board Terra (from 2000) and Aqua (from 2002) — 
with its 36 spectral bands, from 0.47 to 14.24 /im, including 
two new channels 0.405 and 0.550 /im, with a spatial resolu- 
tion of 10 km x 10 km. A detailed description of the instru- 
ments, parameters and data access url is given by Varela et 
al. 2007. Data from ERS-2/GOME and MSG1/SEVIRI have 
not been included on this work because databases are too 
short-term for statistical analysis, but they do offer valuable 
potential as instruments for the future when the database 
increases. The ENVIS AT /SCIAMACHY database has not 
been used in this analysis because it does not provide higher 
spatial resolutions than OMI and MODIS. 

In Appendix I we indicate the official websites to re- 
trieve datasets and the data formats. The correlation anal- 
ysis between AI and AOD with KV is given in Section 6. 

5.3 Using Aura/OMI data 

Here, we show the results of comparing AI provided by 
TOMS and OMI. OMI on the EOS (Earth Observing Sys- 
tems) Aura platform continues the TOMS record for total 
ozone and other atmospheric parameters and can distinguish 
between aerosol types, such as smoke, dust and sulphates. 

Figure \8\ shows the AI values provided by EP/TOMS 
and OMI/ Aura over the Roque de los Muchachos Obser- 
vatory on La Palma. Note the dispersion in the TOMS 
data. This improvement shown in the OMI data (much 
less dispersion) derives from better horizontal and verti- 
cal spatial resolution compared with its predecessor, TOMS 
in Earth Probe. Moreover, TOMS data from 2002 should 
not be used for trend analysis because of calibration errors 
I http:/ /jwocky.gsfc.nasa.gov/nes/news.html). 



5.4 Using Terra/MODIS and Aqua/MODIS data 

MODIS on board the NASA Terra and Aqua satellites pro- 
vides not AI values but an equivalent parameter, the aerosol 
optical depth (AOD). 

In Fig. [9] we show AOD for Terra and Aqua; we see 
that there exists a good relation between the maximum and 
minimum AOD values. The consistency of both data sets is 
excellent. 

The AI threshold for dusty nights is 0.6; we now have to 
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Figure 6. Plumes of dust over the Canary Islands from the west coast of Africa observed by the Moder- 
ate Resolution Imaging Spectroradiometer (MODIS) on NASA's Aqua satellite on 2007 March 10 (full image at 
http://earthobservatory.nasa.gov/Newsroom/NewImages/Images/canary_amo_2007069_lrg.jpg). The easternmost island is just over 100 
kilometres from the African coast. The peaks of La Palma and Tenerife remain clear of dust. 



Table 1. Overview of instruments on board satellites that provide parameters useful for our work 



SATELLITE/INSTRUMENT 


HORIZONTAL RESOLUTION 


PARAMETER 


PERIOD 


Terra/MODIS 


10x10 km 2 


Aerosol Optical Depth (AOD) 


from 2000 


Aqua/MODIS 


10x10 km 2 


Aerosol Optical Depth (AOD) 


from 2002 


Aura/OMI 


From 13x24 km 2 to 24x48 km 2 


Aerosol Index (AI) 


from 2004 



determine the AOD threshold for dusty events, but AI and 
AOD are not in a 1:1 ratio (they depend on refractive in- 
dex, particle size distribution and height of the atmospheric 
layer). 

To determine the AI and AOD threshold for dusty days 
we used information provided in a collaboration between the 
Spanish Environment Ministry, the Upper Council of Scien- 
tific Researches (CSIC) and the National Institute of Mete- 
orology (INM) for the study and analysis of the atmospheric 
pollution produced by airborne aerosols in Spain. They pro- 
vide us with the days in which calima (dust intrusion) oc- 
curred in the Canary islands. As we later demonstrate in 
the plots, in our opinion, these events happened lower alti- 
tudes than those of the observatories so they do not influence 
the measurements of atmospheric extinction. In Fig. 1101 we 
show AOD for Terra and Aqua and AI for OMI. AI values 
are larger than 0.6, and most of AOD data points for Terra 
and Aqua fall above a limit greater than 0.10. This limit is 
consistent with that provided by Romero & Cuevas (2002). 
Therefore, AI>0.6 (0.1 units smaller than provided by Siher 
et al. 2004) and AOD>0.1 will be the thresholds for dusty 
episodes. 



6 RESULTS II: COMPARISON OF AI 

PROVIDED BY OMI AND AOD PROVIDED 
BY MODIS WITH KV FROM THE CAMC 

In this section we analyse the correlation between the aerosol 
index (AI) and the aerosol optical depth (AOD) provided 



by the satellites and the atmospheric extinction coefficient 
(KV) measured by in situ techniques (CAMC). 

6.1 AI from OMI vs KV 

In this subsection we argue about the correlation between 
atmospheric extinction in V (551nm) and the aerosol index 
measured for OMI in the same channels of TOMS (331 and 
360 nm) represented in Figure 1111 In order to determine 
whether there exists a correlation between both parameters 
we also consider the presence or otherwise of clouds and the 
situation for days with calvma. 

In Fig. [TT] we observe an interesting situation not re- 
vealed in the same graph for the aerosol EP/TOMS in- 
dex (Fig. [SJ). Once we put the atmospheric extinction (KV) 
threshold (red line) and the limit for absorbing aerosols 
(green line), the situation reveals four quadrants. In the 
first (the upper right) we have a situation with absorbing 
aerosol and values of atmospheric extinction over the limit 
of photometric dusty nights. In the second quadrant (upper 
left) there are no points that fall inside; this is very impor- 
tant because it means that non-absorbing aerosols do not 
influence the extinction above this threshold. In the third 
quadrant (bottom left) there are still non-absorbing aerosols 
but for clear nights, and this means that this type of sul- 
phate and/or marine aerosol does not influence the extinc- 
tion as we expect. In the last quadrant (bottom right) ab- 
sorbing aerosols are again seen in the presence of low extinc- 
tion values. An explanation could be the presence of weakly 
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Figure 7. Overview of parameters (vertical labels) and periods provided by instruments on board ESA and NASA satellites (shown 
with different colour codes). For the aerosol content only Aura/OMI, Terra/MODIS and Aqua/MODIS results are shown in this paper 
to provide high spatial resolution and a longer-term database. Data from ERS2/GOME and MSG1/SEVIRI have not been included 
in this work because the databases are too short-term for statistical analysis. The ENVISAT/SCIAMACHY database has not been 
used in this analysis because it does not provide higher spatial resolutions than OMI and MODIS. 
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Figure 8. Aerosol Index provided by EP/TOMS (blue points on the left) and Aura/OMI (pink points on the right) in time over the 
Roque de los Muchachos Observatory on La Palma. The improvement shown in the OMI data (much less dispersion than TOMS data) 
derives from better horizontal and vertical spatial resolution compared with its predecessor, TOMS in Earth Probe. 



absorbent particles, such as carbonaceous grains or clouds, 
or a more complex situation with a mixture of cloudy and 
aerosol scenarios, or the presence of absorbing aerosols be- 
low the level of the observatory. The lack of correlation indi- 
cates that AI from OMI should not be used for atmospheric 
extinction characterization at the Roque de los Muchachos 
Observatory. 

We added pink points as dust episodes of calima re- 
trieved at the ground level, yellow being those that corre- 
spond to cloud presence (reflectivity at 331 nm greater than 
15%). This plot is very informative because from the pink 
points we can see that these calima episodes are mainly be- 
low the threshold of dusty nights, this means that they do 
not affect the measurements of extinction because the cal- 



ima is at a lower level with respect to the altitude of the 
observing site, such dusty events reaching the Roque and 
being detected only on some occasions. We also note that 
the clouds are almost all below this limit. 

This means that values above the threshold of dusty 
nights are dust absorbing particles at the height of the ob- 
servatory. We see that in this zone of the plot there are 
almost no clouds and some calima events that reach the as- 
tronomical site. Thanks to strong convective motions, they 
may be driven along the orographic contours of the caldera 
that borders Roque de los Muchachos. In fact, the spatial 
resolution of OMI ranges from 13 km x 24 km to 24 km x 
48 km and can contain part of the caldera. 
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Figure 9. Aerosol optical thickness provided by Terra/MODIS (blue points) and Aqua/MODIS (pink points) in time over the Roque 
de los Muchachos Observatory on La Palma. 
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Figure 10. AI from OMI (blue rhombus) and AOD from MODIS (pink squares correspond to Terra values and green triangles correspond 
to Aqua data) under dusty conditions detected by the INM from 2004 to 2006. 



6.2 AOD from Terra/MODIS vs KV 

The MODIS Aerosol Product monitors the ambient aerosol 
optical depth over the oceans globally and over a portion of 
the continents. The aerosol size distribution is also derived 
over the oceans, and the aerosol type is derived over the con- 
tinents. Therefore, MODIS data can help us to understand 
the physical parameters of aerosols affecting the Canarian 
Observatories. Daily Level 2 data are produced at 10 km 
x 10 km spatial resolution. Aerosols are one of the great- 
est sources of uncertainty in climate modelling. They vary 



in time and in space, and can lead to variations in cloud 
microphysics, which could impact on cloud radiative prop- 
erties and on climate. The MODIS aerosol product is used 
to study aerosol climatology, sources and sinks of specific 
aerosol types (e.g. sulphates and biomass- burning aerosols), 
interaction of aerosols with clouds, and atmospheric cor- 
rections regarding remotely sensed land surface reflectance. 
Above land, dynamic aerosol models will be derived from 
ground-based sky measurements and used in the net re- 
trieval process. Over the ocean, three parameters that de- 
scribe aerosol loading and size distribution will be retrieved. 
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Figure 11. Correlation between AI provided by Aura/OMI and KV from the CAMC. Pink points indicate dust episodes of calima 
retrieved by the INM, and yellow points correspond to cloud presence. Red and green lines indicate the KV and AI limits for dusty nights 
respectively. We can see that most of the points indicating the presence of absorbing aerosols (AI larger than 0.6) but non atmospheric 
extinction (KV smaller than 0.15 mag/airmass) correspond to the presence of clouds and calima below the level of the Observatory. 



There are two necessary pre-assumptions in the inversion 
of MODIS data: the first concerns the structure of the size 
distribution in its entirety and the second one log-normal 
modes that will be needed to describe the volume-size dis- 
tribution: a single mode to describe the accumulation mode 
particles (radius < 0.5 /im) and a single coarse mode to 
describe dust and/or salt particles (radius > 1.0 /im). The 
aerosol parameters we therefore expect to retrieve from the 
aerosol Level 2 product are the ratio between the two modes, 
the spectral optical thickness and the mean particle size. 

The quality control of these products will be based on 
comparison with ground and climatology stations. 

The parameters that we have used from Level 2 in our 
work are: latitude and longitude as geolocation values, and 
AOD (aerosol optical thickness at 0.55 /jm for both Ocean 
[best] and Land [corrected] in a valid range from to 3), 
aerosol type land that contains = mixed, 1 = dust, 2 = 
sulphate, 3 = smoke, 4 = heavy absorbing smoke, and cloud 
fraction land and ocean in percentage. 

In Fig. [12] we show AOD measured by MODIS on board 
Terra against KV. 

In this plot we find two tails of data within a certain 
dispersion of values, and once more the majority of points 
that fall under the threshold of KV smaller than 0.15. The 
flatter tail groups together most of the AOD values lower 
than 0.2, i.e. non-absorbing (< 0.1) or weakly absorbing 
aerosols (marine particles, clouds above ocean and mixed 
scenarios with salt, sulphate particles and clouds). 

In Fig. [12] we distinguish among terrestrial aerosols 
(composed of a mixture of dust and smoke sulphates) and 
marine particles marked respectively by blue and pink points 
and we distinguish the presence of dust events coming from 
Africa (yellow points) using the data provided by INM at 
ground level. 

We see that pink values are situated in the tails of the 
plots and are present mainly at lower KV and AOD values, 
where we expect sulphate sea-salt aerosols and CCN; only 



in some cases are they detected at higher AOD with small 
KV values, perhaps because these aerosols and clouds are 
below the level of the observatory. The second tail covers two 
quadrants corresponding to small AOD but large KV due to 
the presence of weakly absorbent aerosols at the level of the 
Observatory or to a mixture of particles and clouds, and to 
large AOD and large KV, mostly terrestrial aerosols at the 
level of the observatory. The calima (yellow points) is also 
below the KV threshold, perhaps because of their closeness 
to ground level, except for some cases in which it can reach 
the level of the observatory. 

6.3 AOD from AQUA vs KV 

We now perform the same analysis for the MODIS instru- 
ment on board the Aqua satellite that retrieves the same 
aerosol Level 2 data. In the Fig. [13] we show all the Aqua 
data in order to see the correlation between AOD and KV. 
In this case also there is evidence of the two tails with a 
wide dispersion in the data, here, most of the points being 
below the KV threshold. 

We plot terrestrial (blue) and marine (pink) aerosols, 
the yellow points indicating days with calima. The com- 
ments are very similar to those made previously, i.e. the 
marine aerosols are more clustered at small AOD values in 
a lower tail, whereas the terrestrial values are somewhat uni- 
formly distributed with some points at large KV indicating 
absorbing aerosols or clouds at the level of the observatory. 
Most of points fall below the threshold for dusty nights just 
as the majority of calima days. The interesting fact is that 
only the episodes of dust with large KV follow a good linear 
relation with the AOD. 



7 SUMMARY AND OUTLOOK 

From this study we may draw the conclusions listed below. 
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Figure 12. Correlation between AOD provided by Terra/MODIS and the KV measured by the CAMC, distinguishing among the 
terrestrial (mixed dust, smoke, sulphate — blue points), marine (pink points) aerosols and the presence of dust events over the Canary 
islands coming from Africa and collected at ground level (yellow points). The red and green lines indicate the KV and AOD limits 
respectively for dusty nights. 
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Figure 13. Correlation between AOD provided by Aqua/MODIS and the KV measured by the CAMC, distinguishing among the 
terrestrial (mixed dust, smoke, sulphate — blue points), marine (pink points) aerosols and the presence of dust events over the Canary 
islands coming from Africa and collected at ground level (yellow points). The red and green lines indicate the KV and AOD limits 
respectively for dusty nights. 



We have compared the AI measurements provided by 
the TOMS on board Earth Probe satellite (Level 3) with the 
atmospheric extinction coefficient provided by the CAMC at 
the ORM (2400 m above MSL). The main causes of the lack 
of correlation between both parameters are: 

• The TOMS Level 3 data considered in this paper have 
a resolution of 1° x 1.25° , so the AI is averaged over areas 
whose size covers the entire islands of La Palma and Tener- 



ife. High resolution Level 2 data should be used for a better 
fit (IFOV of 35 km x 35 km). 

• The TOMS is very sensitive to the presence of highly 
reflective clouds because it uses channels centred on the UV 
to measure AI. Moreover, AI incorporates absorbing parti- 
cles in ranges that do not affect atmospheric transparency 
in the visible range. 

• TOMS measurements are retrieved at local noon while 
CAMC values are averaged over night hours. 



Astronomical site selection: On the use of satellite data for aerosol content monitoring 13 



For this reason the EP/TOMS database is not useful for 
the characterization of the presence of dust above either the 
Canarian astronomical observatories (2400 m above mean 
sea level) or for other high mountain sites (Mauna Kea and 
San Pedro Martir). 

We have explored the use of other detectors on board 
different satellites that operate in bands of astronomical in- 
terest (the visible and NIR) and with better spatial resolu- 
tion than TOMS. The selected parameters were the aerosol 
index provided by Aura/OMI — with visible and ultraviolet 
channels and with a spatial resolution from 13 km x 24 km 
to 24 km x 48 km — and the aerosol optical depth (AOD) 
provided by Terra (from 2000) and Aqua (from 2002) in 
MODIS— with its 36 spectral bands, from 0.47 to 14.24 fim, 
including two new channels at 0.405 and 0.550 /im, with a 
spatial resolution of 10 km x 10 km. In order ot obtain the 
best spatial, spectral, radiometric and temporal resolutions, 
we have decided to work only with Level 2 data that have 
the same resolution as the IFOV satellite. 

We conclude that the OMI instrument detects aerosol 
presence with more precision than TOMS and does not de- 
tect non-absorbing particles with high atmospheric extinc- 
tion values (larger than 0.15 mag/airmass). This fact coin- 
cides with expectations because non-absorbing aerosols such 
as sulphates or marine aerosols do not give high extinction 
values (threshold greater than 0.15 mag/airmass). We can 
see that most of points fall at lower extinction values below 
the threshold for dusty nights, suggesting presence of non- 
absorbing or weakly absorbing (e.g. carbonaceous) aerosols. 

In order to obtain the limits for dusty episodes on 
the AOD scale, we have checked the calima days from the 
records of the Instituto Nacional de Meteorologfa de Ca- 
narias Occidental, following NAAPS, ICOd/DREAM or SK- 
IRON models, and we obtain a threshold near AOD>0.1 
units and AI>0.6 for dusty episodes. 

We study where the calima events and cloud presence 
fall in the plot of correlation between atmospheric extinc- 
tion and AI. Dust episodes measured at ground level are 
mainly below the threshold for dusty nights on the atmo- 
spheric extinction scale (KV<0.15 mag/airmass), meaning 
that the presence of calima affects low altitudes, and that 
only in a few cases does it reach the Roque. We also see that 
all clouds detected for their high reflectivity (greater than 
15%) are below the threshold for dusty nights; in only two 
cases are they above 0.15 unit in mag/airmass and in such 
cases they do not correspond to calima events. 

Study of Terra/MODIS data has shown that a great 
number of points fall in a range below 0.40 units for AOD 
and 0.15 mag/airmass for KV, and that two tails are evident: 
the first one has high AOD values for low KV and the second 
has high AOD and high KV, showing a large linear correla- 
tion among both parameters. Is important to underline that 
there are only a few points with low AOD values and high 
KV (data are decontaminated of cloud presence). Chlorides 
and marine aerosols can be well identified and normally do 
not affect the KV and correspond to AOD < 0.2, so this 
means that marine and sulphate aerosols are not absorbent 
as we expect. We can see that the most of points corre- 
sponding to dust (calima) events fall below the threshold of 
0.15 mag/airmass because they are detected near the sur- 
face; only in certain cases can they reach the level of the ob- 
servatory, carried by wind or convective motions, providing 



AOD larger than 0.1 units and KV larger than 0.15. These 
measurements correspond mostly to terrestrial aerosols. 

The study of the Aqua data produces results almost 
identical to those of Terra. In this case there are also values 
at high AOD and low KV and values at low AOD and KV. 
This wide dispersion explains the lack of correlation between 
both parameters. The most populated tail falls below the 
threshold of 0.15 unit in mag/airmass. 

Marine aerosols are more clustered at low AOD values 
(< 0.20) and low KV (<0.15), whereas terrestrial aerosols 
fall in the AOD>0.2 zone. Dusty (calima) days correspond 
mostly to the presence of terrestrial aerosols and are present 
near ground level (inside the area of 10 km x 10 km), so they 
appear at low KV. Only in some cases do they reach higher 
altitudes and become detectable from the astronomical ob- 
servations (large KV); this is the only case that presents 
a linear correlation between both parameters. We therefore 
also need in situ data to distinguish between both situations. 
We must explore other clues for the vertical aerosol drainage 
analysis (winds, humidity, etc.). 

At present, the AI and AOD values provided by the 
NASA satellites are not useful for aerosol site character- 
ization, and in situ data are required to study drainage 
behaviour, in particular at those astronomical sites with 
abrupt orography (ORM, Mauna Kea or San Pedro Martir). 
Spatial resolution of the order of the observatory area will 
be required in these cases. 

Moreover, in order to obain much better spatial reso- 
lution we are now exploring the use of SEVIRI-MSG2 (1.4 
km x 1.4 km) (December 2005-12) for Europe and Africa, 
and in the future, ATLID(LIDAR)-EARTHPROBE (with a 
horizontal sampling interval smaller than 100 m) (2012-15) 
for global coverage. The CALIPSO satellite will be used for 
measuring the vertical structure (drainage) and properties 
of aerosols (the horizontal and vertical spatial resolution are 
5 km and 60 m respectively). 

Ground measurements will be complemented by LI- 
DAR data (INTA) of 30 m resolution (NASA MPL-NET- 
AERONET) and by the IAC airborne particle counter (from 
Pacific Scientific Instruments) installed at the ORM in 
February 2007 (with six channels: 0.3, 0.5, 1, 3, 5, 10 /xm) 
and by the INM Multifilter Rotating Shadowband Radiome- 
ter (MFRSR) programmed to be installed at the ORM in 
the near future (consisting of six narrow passbands between 
414 nm and 936 nm) that will provide the size, density and 
vertical distribution of the aerosols. 

Tests with AERONET in situ data during daytime 
could be very interesting to make a comparison with remote 
sensing satellite data. 
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APPENDIX I 

FORMAT OF SATELLITE DATA 

In this appendix we describe the format of the data for the 
different satellites used in this work and we indicate the urls 
from which to retrieve the datasets. 



1.1 N7/TOMS (1984-1993) 



Level 



data 



available 



at 



http://daac.gsfc.nasa.gov/data/dataset/TOMS/Level_2/N7/ 



and were ordered with ftp-push. NASA loaded the data 
onto our ftp url. These are raw data in hdf5 format that 
must be seen with HDFview and afterwards processed with 
computer programs. Level 3 aerosol data are available at 
ftp://toms.gsfc.nasa.gov/pub/nimbus7/data /aerosori and 
are already available in ftp format. These are daily average 
data with a resolution of 1.25° in longitude and 1° in 
latitude. These data are in ASCII format with 288 bins 
in longitude, centred on 179.375 W to 179.35 E, every bin 
has a 1.25° step. In latitude there are 180 bins centred on 
89.5 S to 89.5 N with a step of 1°. The values are in 3 
digit groups, missing data being tagged 999 and the other 
numbers being multiplied by 10. 

1.2 EP/TOMS (1996-2005) 

Level 2 data are available at 

http://daac.gsfc.nasa.gOv/data/dataset/TOMS/LeveL2/EP/l 
The format and the data processing is the same 
as for Nimbus7. Level 2 data are produced with 
a spatial resolution of 39 km x 39 km at nadir. 
Level 3 aerosol data are immediately available at 
ftp: / / toms.gsfc.nasa.gov / pub /eptoms / data/aerosol . 

It is important to underline that the aerosol monthly 
average datasets are computed using only positive values 
(i.e. absorbing aerosol indices) of the aerosol index for each 
month. Values of zero are used in the averaging whenever the 
aerosol index is negative. The final monthly average datasets 
contain aerosol index values greater than or equal to 0.7. 



1.3 Aura/OMI (2004-NOW) 



Level 
13 x 



data with 
24 km at 



a spatial 
nadir are 



resolution 
available 



of 
at 



|http: //daac.gsfc.nasa.gov/data/dataset/OMI/Level2/OMT03/ 
They can be ordered via ftp-push in the same 
way as the others. To obtain information about 
effective cloud pressure and fraction data go to 



aerosol optical thickness, as well as ancillary information 
produced from the TOMS Version 8 algorithm applied to 
OMI global mode measurements. In the global mode each 
file contains a single orbit of data covering a width of 2600 
km. Compared to TOMS, OMI's smaller field of view results 
in a larger "sea glint" per unit field of view and a correspond- 
ingly larger error in derived ozone under these conditions. 
The OMT03 aerosol index is not valid for solar zenith angles 
greater than 60°. Because the OMI solar zenith angles are 
typically higher than the solar zenith angles for TOMS at 
the same latitude, the OMI AI becomes invalid at somewhat 
lower latitudes than TOMS. This may show a cross-track de- 
pendence in the OMI AI and is not corrected by the radiance 
measurement adjustments (error in the AI up to 4%) . Com- 
pared to TOMS, the OMT03 Aerosol Index is 0.5 NVALUE 
high. Users of the aerosol index are advised to make this 
correction for consistency with the TOMS data record. For 
users not interested in the detailed information provided in 
OMT03 dataset several gridded products are being devel- 
oped. Initially, DAAC will grid OMT03 data in a format 
identical to that used for TOMS (1° x 1.25° lat/long) and 
will make it available through the TOMS website. However, 
to take advantage of the higher spatial resolution of the OMI 
products DAAC intends to produce higher resolution grid- 
ded products for all OMI datasets, including OMT03. 



1.4 Terra/MODIS (2000-NOW) 



data 



are 



available 



at 



Level 2 

http:/ /daac.gsfc.nasa.gov/data/dataset/MODIS/02_Atmosphere/01_Level_2 i 



We can order them with ftp-push. These are raw data in 
hdf5 format and so we must follow a similar procedure to 
the one we use with TOMS. Daily Level 2 (MOD04-Terra) 
aerosol data are produced at the spatial resolution of 
10 x 10 km at nadir. We can also retrieve data via ftp 
(they are already available on the web) through the url 



|ftp://g0dps01u.ecs.nasa.gov/ MODIS_t erra-Atmosphere/M OD04_L2.004 

Level 3 data are available at 
http: / / disc.sci.gsfc.nasa.gov / data/dataset /MODIS /02_ Atmosphere/02_Level 
These are different atmospheric data, daily, weekly and 
monthly averaged in a global 1° x 1° grid. The method to 
order this dataset is ftp-push: The aerosol information is 
stored in MOD08_D3, MOD08_E3 and MOD08_M3. Note 
that in MODIS there is no aerosol index but only the 
aerosol optical thickness. 



1.5 Aqua/MODIS (2002-NOW) 

The same as for Terra/MODIS. 



The urls are: 



http:/ /daac. gsfc.nasa.gov/data/dataset/MODIS-Aqua/02_Atmosphere/01_L 
(the ftp address ftp:/ /gOdpsOlu. ecs.nasa.gov/MODIS_Aqua_Atmosphere/MY 
for Level 2 data and http:/ /disc. sci.gsfc.nasa.gov/data/dataset/MODIS-Aqua 
for the Level 3 data. 
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AEMET Agenda Estatal de Meteorologi'a 
AEROCE Atmosphere/Ocean Chemistry Experiment 
AERONET Aerosol Robotic NET work 



http://daac.gsfc.nasa.gov/data/dataset/OMI/Level2/OMCLDRRXl Aerosol Index 

OMT03 provides aerosol index, total column ozone and AOD Aerosol Optical Depth 
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AQUA Earth Observing System Post Meridian (PM) 
ASCII American Standard Code for Information Inter- 
change 

ATLID ATmospheric LIDAR 

AURA Earth Observing System Chemistry mission 
CALIPSO Cloud-Aerosol Lidar and Infrared Pathfinder 
Satellite 

CAMC Carlsberg Automatic Meridian Circle Telescope 

CCN Cloud Condensation Nuclei 

CSIC Consejo Superior de Investigaciones Cientfficas 

ELT Extremely Large Telescopes 

ENVISAT ENVIronmental SATellite 

EP Earth Probe 

ERS European Remote Sensing Satellites 

ESA European Space Agency 

EUMETSAT EUropean METeo SATellite 

ENVISAT ENVironmental SATellite 

FP6 Sixth Framework Programme 

FTP File Transfer Protocol 

GMT Greenwich Mean Time 

GOME Global Ozone Monitoring Experiment 

GSFC Goddard Space Flight Center 

HDF Hierarchical Data Format 

HTTP Hyper Text Transfer Protocol 

IAC Instituto de Astroffsica de Canarias 

IAU International Astronomical Union 

ICoD/DREAM Dust Loading model forecast from Insular 

Coastal Dynamics 

IFOV Istantaneous Field of View 

INM Instituto Nacional de Meteorologfa 

INTA Instituto Nacional de Tecnica Aerospacial 

IP Internet Protocol 

L2 Level 2 data 

L3 Level 3 data 

KV Atmospheric extinction coefficient in V-band 
LIDAR Light Detection And Ranging 
MFRSR MultiFilter Rotating Shadowband Radiometer 
MML Maritime Mixing Layer 

MODIS Moderate Resolution Imaging Spectroradiometer 
MPLNET MicroPulse Lidar NETwork 
MSG Meteosat Second Generation 
MSL Mean Sea Level 
N7 Nimbus-7 

NAAPS Navy Aerosol Analysis and Prediction System 
NASA National Aeronautics and Space Administration 
NEODC NERC Earth Observation Data Center 
NERC Natural Environment Research Council 
NIR Near Infra Red 
NRT Near Real Time 
NUV Near Ultra Violet 
03 Ozone 

OAI Izana Atmospheric Observatory 
OMI Ozone Monitoring Instrument 
0MT03 OMI Total Column Ozone 

OPTICON Optical Infrared Coordination Network for 
Astronomy 

ORM Roque de Los Muchachos Observatory 
OT Teide Observatory 

SCIAMACHY Scanning Imaging Absorption SpectroMeter 
Atmospheric Chartography 

SEVIRI Spinning Enhanced Visible and InfraRed Imager 
SKIRON Weather Forecasting Model operated by Univer- 



sity of Athens 

TERRA Earth Observing System Anti Meridian (AM) 

TL Troposphere Layer 

TOMS Total Ozone Mapping Spectrometer 

UK United Kingdom 

UV Ultra Violet 
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